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Abstract 

From microscopic models, a Langevin equation can in general be derived only as 
an approximation. Two possible conditions to validate this approximation are 
studied. One is, for a linciar Langevin equation, that the frequency of the Fourier 
transform should be close to the natural frequency of the system. The other is 
by the assumption of 'slow' variables. We test this method by comparison with 
an exactly soluble model, and point out its limitations. We base our discussion 
on two approaches. The first is a direct, elementary treatment of Senitzky. The 
second is via a generalized Langevin equation as an intermediate step. 
PACS: 05.10.Gg, 05.30.-d, 05.40.Ca 

1 Introduction 

The Langevin equation is certainly not exact (except in the special case when 
the underlying equations of motion are all linear), and so it is of interest to find 
the conditions for it to be a good approximation. 

The equation has been justified, starting from microscopic equations of mo- 
tion, by various general arguments and in some particular models. These argu- 
ments start from the Heisenberg equations of motion for a simple (macroscopic 
or mesoscopic) systcmi in interaction with an environment (or bath) with many 
degrees of freedom. We will assume the systems are quantum ones, but every- 
thing we say would equally apply to classical systems (replacing commutators 



by Poisson brackets, etc.) We choose a microscopic model in which the interac- 
tion is bihnear in the environment variables. This is the simplest change form 
the (trivial) linear case, and is necessary if the environment is fermionic. 

We will concentrate on two approaches. The first is due to Senitzky [I]. It 
uses only elementary quantum mechanics, but it is not clear for what ranges of 
parameters the approximations made are good ones. . We argue that one region 
in which Senitzky's approximations may be justified is where the frequency is 
close to the natural frequency of the simple system (which therefore assumes 
that system to be linear). Our approach is to regard the Langevin equation as 
summing an infinite subset of perturbation theory terms. The question is then, 
when does this subset dominate? 

For an earlier discussion of Senitzky's argument see [TT]. So far as we know, 
text books do not offer a derivation as simple and direct as Senitzky's. See for 
example [3], g], [5], [6], [7], 0. 

The second method of justification is indirect. First a generalized Langevin 
equation is established by a projection method (TH], [1]. This is formally 
exact, but probably not useful. (It is usually restricted to the case of a linear 
system, equation (2.3).) In section (5.2) we point out that in the generalized 
Langevin equation the separation into a noise term and a dissipative term is 
ambiguous. 

Then it is assumed that the system variables are 'slow' compared to the 
environment ones. It is argued that the generalized Langevin equation then 
simplifies greatly, and becomes an ordinary Langevin equation. In section 5, we 
examine this proposal critically, and test it in an exactly soluble example. For 
some text book accounts, see [8], [2], |3- 

It is worth noting that both the above methods work by manipulating the 
underlying Heisenberg equation of motion. Yet the Langevin equation in the 
end refers to an expectation value over some distribution function, chosen on 
physical grounds. Clearly some assumptions have come in during the course of 
the derivations. 

2 The underlying microscopic systems 

We use the Heisenberg picture throughout. 

We are concerned with attempts to deduce a Langevin equation from an un- 
derlying dynamics. This consists of a simple macroscopic or mesoscopic system 
interacting with an environment (sometimes called a bath). For simplicity, we 
take the system to have one degree of freedom, with phase space Q, P. The 
environment has many (microscopic) degrees of freedom qi,Pi with i — 1, N. 
The Hamiltonian is 

H = Hoit)+H'{t), Ho{t)^Hs{t)+HEit) (2.1) 

N 

Hs = {l/2)P^/M+V{Q). He - Y.^\/2){pl /m,+m,ujUl), H' = -aQK{q„p,) 

i=i 

(2.2) 



2 



{N large). Thus we are assuming the environment to be a set of oscillators. We 
will sometimes take 

V{Q) = (M/2)02q2 (2.3) 

that is the system is an oscillator too. In this case, the resulting Langevin equa- 
tion is linear. In (2.2), a. designates the coupling strength. We have assumed 
that H' is linear in the system variables, and for simplicity chosen Q (not P) 
to appear there. 

It will be useful to define operators Qq, Pq, qai,pai to coincide with Q, P,qi,pi 
at an initial time to, but to vary with time according to the free {a = 0) equations 
of motion. That is 
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Qo(to) = Q{to)- Po{to) = P{to), 
hQo{t) = i[Hso, Qo{t)], hPo{t) = i[Hso, Po{t)], 
Hso = Hs{Qo{t),Po{t)), 

qOi{to) = Qiito), Poi{to) =Pi{tQ), 

Moi{t) = i[HEo, qoz{t)], fipoiit) = i[HEO,Poi{t)], 

HEO = HE{qo^{t),Po,{t)). (2.4) 

Note that Hso and Heo are each independent of t, 

If the function JsT in (2.2) is linear, the equations of motion of the environment 
are hnear, and these variables can be simply eliminated to produce an exact 
Langevin equation. No statistical distribution function appears in this Langevin 
equation. For these two reasons, the case of K being linear is misleadingly 
simple, although it is treated in many text books, for example as a model of 
Brownian motion and of decoherence. 

The simplest non-trivial example is for K to be bilinear in the qi,Pi, and 
this is what we shall assume below. Note that the environment might consist 
of fermionic variables (like conduction electrons), and then K would necessarily 
be bilinear or of higher degree. We will not treat this fermionic case explicitly, 
but our arguments below can easily by generalized to cover fermi statistics. 

The remaining ingredient is a statistical distribution function. In non- 
equlibrium statistical physics, the choice of distribution function is a matter 
of physical judgement. We shall mention just two possibilities. The first is 
simply 

p = Z-^ exp{-H/T) (2.5) 
the equilibrium distribution. The second is the factorized free distribution 

po = Z^^ cxp{-Ho/T) ^psxpE^ PS X Ze' cxp{~Heo/T) (2.6) 

where Hso is the time-independent energy defined in (2.4), and we need not 
specify pE further. Units of temperature T are chosen so that Boltzmann's 
constant is unity, and the partition functions Z, Zq and Ze are normalization 
factors. The use of po is motivated by the idea that the environment is initially 
at equilibrium by itself, and then the system is brought into contact with it at 
some initial time to. Since we use the Heisenberg picture, density matrices are 
time-indpendent, and so the factorization property (2.6) is a single condition, 
not one for each value of t as it would be in the Schrodinger picture. But 
Heo, and therefore also pE, depend implicitly on the initial time to. For the 
Hamiltonian He in (2.2), pEo factorizes further: 

pEO = Y[pio (2.7) 

i 

We define expectation values, for an operator X, 

{X)=tT{Xp), (X)o = tr(Xpo). (2.8) 
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In its simplest form, the Langevin equation which one might hope to derive, 
is usually assumed to be (for t > to) 

POO 

MQ{t) + V'{Q{t)) + aM dt'Cit - t')Q{t') = aKo{t), (2.9) 

J to 

where 

Cit, t') = {C{t, t'))o = {C{t, t'))E = tv{CE{t, t')pE} (2.10) 

and 

C{t,t') = e{t-t')c{t,t'), c{t,t') = -^[Ko{t),Ko{t')] = C{t,t') - C{f'A) 

(2.11) 

and K{) = K{qoi,pai). (In the classical case, the right hand side of (2.11) is 
to be replaced by the Poisson bracket and the trace in (2.8) by J fl^dqidpi.) 
Equation (2.9) shows the characteristic features of a Langevin equation: the 
noise Kq and dissipation contained in the C term (non-Markovian in general, 
that is frequency dependent.) 

If K is linear in the environment variables, (2.9) is exact .This follows simply 
by solving for the qi{t) (the retarded solution) in terms of Q and qi{to),pi(to) and 
inserting this solution into K(qi,pi). In this special case C{t,t') is a c-number, 
and in fact is not dependent of any dynamical variables; so the expectation 
value in (2.10) is redundant. 

A typical quantity one might want to compute by using (2.9) is the correla- 
tion function 

Sqit, t') = {Q{t)Q{t'))o = tT{Q{t)Qit')po}. (2.12) 

In general, Sq may not be a function of {t — t') only. This is because po defined 
in (2.7) does not commute with H' and so not with the total Hamiltonian H. 

If t and t' are much later than to, the form of ps may be unimportant in (2.12). 
There is a general result, the fluctuation-dissipation theorem. Define 

SKoit - t') = {Ko{t)Ko{t'))E = tv{Ko{t)Ko{t')pE} (2.13) 

(note that 5*^ does depend only on the time difference because UKq = i[HE, Kq] 
and [He,Pe] = 0) and define the Fourier transform by 

Sko = j doje-"^^*-''^SK{oj). (2.14) 

Then the relation is 

-inc{u) = tanh (^^^ (5k(w) + SK{-io)) (2.15) 

where c is the Fourier transform of the expectation vlue of c in (1.8). 

In section 3 we review Senitzky's derivation of a Langevin equation, and then 
derive one condition for its validity. Section 4 reviews the so-called generalized 
Langevin equation, and section 5 discusses how the ordinary Langevn equation 
might follow if Q, P are 'slow' variables. Section 6 summarizes our conclusions. 
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3 Senitzky's argument 

We first emphasize the sahent features of Senitzky's [T] argument, which is quite 
general. 

By using the Heisenberg equations of motion, for K and He^ Senitzky derives 
the exact equations 

MQ + V'{Q) = aK, (3.1) 
K{t)=Ko{t)-— / dh / dt2U*{t-ti)[K{h),K{t2)Q{t2)]U{t-h), (3.2) 

" Jto Jto 

where 

ihkit) = [K{t),H] = [K{t),HE{t) + Hs{t) + H'{t)] = [K{t),HE{t)], (3.3) 

U{t) = c^j>{-itHEito)/fi) = exp{-itHEa/h). (3.4) 

(Note that it is HE{t) which comes in (3.3), not Heq ) 

Senitzky then approximates in (3.2) K by Kq and U*{t — ti)Q{t2)U{t — ti) 
by Q{t2), to obtain 

ia /■* Z"*^ 

K{t)^Ko{t)-— dh dt2[K„{t),Ko{t2+t~h)]Q{t2). (3.5) 

Jto Jto 

Now the tl integration may be done to give 

Kit) Koit) - a [ dt'C{t,t')Q{t'), (3.6) 

Jto 

where C is given by (2.11). 

Equation (3.6) is the first stage of Senitzky's approximation. The second is 
to use to use the further approximation 

C{t,t')^C{t,t'), (3.7) 

giving, together with (3.1), the Langevin equation (2.6). (Note that the approx- 
imation (3.7) is needed in the classical case as well as the quantum one. Only 
in the linear case is (2.9) exact.) 

The approximations leading to (3.6) amounts to neglecting in K some, but 
not all, terms down by a power of a^. In order to discuss this, we define the 
power series 

K = Ko + aKi+a^K2 + Q = Qo + aQi + a^Q2 + ■■■■ (3.8). 

Then for example a neglected term containing KiQq is not obviously smaller 
than the retained term KqQi in (3.5). Thus it is far from obvious that (3.5) is 
a valid approximation. 
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3.1 Perturbation theory 



In order to investigate the region of validity of Senitzky's first approximation 
(3.5), we consider an expansion of K in powers of a. The approximation picks 
out an infinite subset of terms. The question is, when do these terms dominate. 
Wc shall attempt to answer this question by looking at the lowest order terms, 
but we believe that our argument generalizes to all orders. 

The solution of the exact equation (3.4) is more easily derived direcly from 
the original Heisenberg equations of motion (3.3) etc. It is 

K{t)=W*{t)Ko{t)W{t), (3.1.1) 

where W is the solution of 

iTiWit) = H'o{t)W{t) = -aKo{t)Qoit)W{t), W{to) = 1. (3.1.2) 

We may now compare the order term, K2, from (3.1.1) with the approx- 
imation from (3.5). To order a^, the exact (3.1.1) gives 

^2(t) = -4/ dh f dt2[[Ko{t),Ko{h)Qoiti)],Ko{t2)Qo{t2)]. (3.1.3) 

Ti Jta Jto 

In order to find the approximate form of K2 deduced from (3.5), we first 
find Qi . This is derived (3.1), with K on the right approximated by Kq. The 
solution of this equation may be written 

Qi{t) = I ^* dt'[Qo{t),Qo{t')]Ko{t'). (3.1.4) 

(This is easy to check if V is a sum of integral powers of Q.) 
Inserting (3.1.4) into (3.6) gives (using (2.11)) 

K2{t) ^ -i / dh f dt2[Koit),Ko{ti)]Ko{t2)[Qo{h),Qo{t2)]- (3.1.5) 

il' Jto Jto 

The difference between (3.1.3) and (3.1.5) is 

-i / dh [ ' dt2[[Ko{t),Ko{h)],Ko{t2)]Qo{t2)Qo{ti). (3.1.6) 
n Jto Jto 

In general, (3.1.5) and (3.1.6) are of the same order. It is clear that the 

Langevin equation sums a subset of the terms in (3.1.1). The question is: what 
are the conditions for this subset to dominate over the other terms, like (3.1.6)? 
We are able to propose an answer to this question, but only for the case that 
the system is an oscillator, that is assuming (2.3). Then the free equation for 
Q{t) has a Greens function 

G{t) = e{t){Mn)-^ sin(m). (3.1.7) 
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and 

0{t - t')[Qo{t), Qo{t')] = -ihG{t - t'). (3.1.8) 
The Fourier transform of (3.1.7) is 

The structure of (3.1.5), but not of (3.1.6), allows for a pole in (3.1.9) at w = f2 
to appear in suitable functions of co. But to complete this argument, we must 
include Senitzky's second approximation, and to do this we need a more detailed 
model for K, which is the subject of the next subsection. 



3.2 A model interaction 

In this subsection we will take to = — oo. We discuss later what changes if to is 
finite. 

In order to test Senitzk's argument, we take the simplest model we can 
think of, which is not totally linear. We take K in (2.2) to be bilinear in 
the environment variables. If we have in mind a model for Brownian motion, 
then the simple case when K is linear might be thought of as representing the 
Brownian particle emitting phonons. In our bilinear generalization, phonons 
scatter off the particle. Of course, these are only crude models. 

Our formalism can easily be extended to a fermionic environment, when K 
cannot be linear but can be bilinear. 

Define annihilation operators by 

ipi{t) + miU>iqi{t) f ■ s\ 

a^{t) ^ , aj(i) = exp(-«w^i)aj. (3.2.1) 

Then we take 

N 

K = k-{k)E where k{t) = L"^ ^ {u:,ujjf'^a*{t)bijaj{t), (3.2.2) 

where the bij are dimensionless numbers and i is a parameter with dimensions of 
length which we shall not specify further (but which might represent for example 
a lattice spacing). With K defined as in (3.2.2), a in (2.2) is dimensionless. 
We shall concentrate on c, the odd part of C in (2.11). Prom (3.2.2), 

c{t,t') = '^{uJiOjjy/^Xij{t,t')a*aj (3.2.3) 

where 

Xij{t,t') = -iL~^^^uikbikbkj[exp{i{uii-uJk)t+i{(jJk-oJj)t'}-{i ^ j)]. (3.2.4) 
k 
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From (3.2.3), (3.2.4), 



c{t - t') = {c{t,t'))E = J2^iM^i) 



(3.2.5) 



where 



Xii{t - t') = -2L ^ ^ ojkhkhi sin[(w/c - uji){t - t')] 



(3.2.6) 



k 



and 



w{uj) = fi.w[exp(?lw/T) — 1] 



-1 



(3.2.7) 



The classical (high temperature) limit of w{uj) is T, and the low temperature 
limit is huexp{—?iu/T). 

Wo first discuss the Scnitzky's second approximation (3.7), and return to 
(3.5) later. So we begin by assuming (3.6) and combine it with (3.1) to give 



(Wc have symmetrized the order of the operators in the integrand in order to 
make it explicitly Hermtian. Since [C, Q] = 0{a) this involves only higher order 
terms in the integrand.) 

We will formally solve (3.2.8) as a power series in a, and then compare the 
terms in this series, with and without the use of (3.7). Of course, the Langevin 
equation sums an infinite number of orders of a. but approximations which fail 
term-by-term in the power series are unlikely to succeed for the complete series. 

In order to be definite, we assume we are using the Langevin equation in 
order to calculate the Q-noise correlation function (2.12). It is this assumption 
which gives a physical motivation for using po as the distribution function. Up 
to this point, there was no reason to prefer pq in the approximation (3.7). 

Through order a^, we get 



Sgit - t') = {Qo{t)Qo{t'))o + a' / dhdt2G{t - h)G{t' - t2){Ko{h)Kc,{h))o 



(Note that, using tg = — oo, Sq turns out to depend only on t — t', athough this 
is not obvious from the definition (2.12).) If the approximation (2.7) were valid, 
we should be able to approximate the expectation value at the end of (3.2.9) by 




(3.2.8) 




x{{C{tut2)Koit3) + Ko{h)C{h,t2)}Ko{U))o 
+(herm. conj. witht -H- 1') 



(3.2.9) 



(3.2.10) 



etc. So this is what we now check, using (3.2.3) and (3.2.4). 
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To evaluate the expectation value in (3.2.10), we need (using (3.2.2)) 



J{prA. (3.2.11) 

ijklmn \ t / 

This receives nonzero contributions from the following values 

{&) : i = j, k = n, I = m; {h) : i = I, k = n, j = m; {c) : i = n, j = k, I = m. 

(3.2.12) 

where in each case we assume there are no other equahties. There arc also 
contributions like i=j = k = n, l = m, but for large these are negligible 
compared to (3.2.12). 

The contribution to (3.2.11) from region (a) in (3.2.12) gives the right hand 
side of (3.2.10); so the question is, when are the contributions from (b) and (c) 
small compared to that from (a)? 

The contribution to C from (3.2.9) has the form 

ia^[G{u})]'^ Y^[Wa + Wb + We] + herm.conj. (3.2.13) 

ijkl 

where Wa, Wb, Wc come from the ranges in (3.2.12): 

Wa = G{uj)w{ujj)w{uji)bkibikhjbji{{uj - + Uj + ze)"M(w - Wfe + wj) 

-{^i,j,k,i -<^i,j,k,i)^ (3.2.14a) 
Wb = G(uJi - u)k)w{iOj)w{u)i)bijbjibikbki{LO - + wj + ie)~^(5(w - + cjj), 

-{oJi,j,k,i -i^i,j,k,i), (3.2.146) 
Wc = G{uJj - uJi)w{ujj)w{uJi)bijbjibikbki[{}^ - + ujj + ie)~^5{uj - Wfc + loi)) 

-{(^i,j,k,l -^^i,j,k,l), (3.2.14c) 

Here G is defined in (3.1.9), and has a pole where its argument is equal to 
the natural frequency .Tthere are three such poles at w = Q in (3.2.13) 
from(3.2.14a), whereas from the other two terms there are only two poles. Thus 
Wa, and so the approximation (3.2.10) to (3.2.9), may be good for values of uj 
sufficiently near to Q. Otherwise, we can see no reason why Wa, Wb, Wc should 
not be comparable. 

The structure of the terms in equations (3.2.14) is illustrated in Fig. 1. The 
generalization to all orders of a is that the Langevin equation sums 'bubble 
graphs', that is graphs like the first in the figure, with a sequence of 'bubbles' 
like the one marked X in the figure, connected by Q-propagators. These graphs 
have the maximum number of poles at w = il. 

We now turn to the first approximations made in Senitzky's argument, the 
replacement of (3.2) by (3.5). To order a^, this implied the neglect of (3.1.6) 
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Figure 1: Graphs symbolizing the teri):ig in (3.2.14), Thick Unes represent the 
Q-propagator (3.1.9). Thin lines indicate the pairings in (3.2.12). The X repre- 
sents the commutator c in (2.11). and black circles represent the Kq operators in 
(3.2.9). The vertices are labeled according to the values of the times in (3.2.9). 



compared to (3.1.5). If we work out the contribution of (3.1.5) to (3.2.9), we 
need instead of (3.2.10): 



{[C{h,t2),Ko{ts)]Ko{t4))o. 



(3.2.15) 



A contribution like (3.2.14a) it is zero, because [C{ti,t2), KQ{t^)] = 0. Thus 
(3.1.5) has parts which with three poles at u; = f2, but (3.1.6) does not, so near 
this pole both Senitzky's approximations (3.5) and (3.7) are justified. 

If to is finite, there are no poles at lu = fl. However, in the Fourier transform 
of (3.2.9), times less than to do not contribute much if 



Thus we expect the would-be poles at a; = O to be large provided that |<:o| is 
large enough to satisfy (3.2.16). 

4 Non-linear Q-systems 

Most of our arguments above have assumed that the simple system is an oscilla- 
tor, as in (2.3), with a well-defined natural frequency fl, so that we can recognize 
the poles (3.1.9). But we may ask if the arguments extend to the more general 
potential V{Q) in (2.2), when the equation of motion for Q is nonlinear. Then 
a simple Greens function like (3.1.9) does not exist. However, when working to 
first order, there is the equation (3.1.4), where 0{t — t')[QQ{t), Qo(i')] looks a bit 
like a Greens function, though in general it is an operator and is a function of 
t and t' separately. 

Let us consider an example like (3.2.9). For a non-linear system, each Greens' 
function G{t — t') is replaced by the operator i9{t — t')[Qo{t),Qo{t')]/h. We 
have to take the expectation value of products of these commutators using ps 
in (2.6). In order to make the same sort of argument as we did with the prop- 
agator poles (3.1.9), two things are needed: the expectation values of products 
commutators should approximately factorize into products of expectation values 
(Qo(i)Qo(^')); '^.nd the Fourier transforms of each of these should have a pole 
at some definite frequency. Whether or not these things happen in any useful 
approximation will depend upon the potential V{Q). 

As the simplest example, take a product of two commutators, as would occur 
in the first line of (3.2.9). In the nonlinear case, the product of two Greens 
functions G{t — t{)G{t' — t2) would be replaced by 



-{i/n^)eit - h)e{f' - t2M[Qoit), Oo(ii)][0o(t'), Qoit2)]ps}, (4.i) 



Let us assume that Hs has a discrete set of energy eigen-states \a) with energy 
Ea- and take the simplest case where ps corresponds to the (pure) ground state 
|0). Then (4.1) has contributions 



-{l/h')e{t-t^)eit'-t2) J2 ((0|Qo(i)|a)(a|Qo(ti)|^)(/3|Qo(i')l7)(7|Qo(t2)|0) 



U)\to\ > 1 



(3.2.16). 
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-(< o h) - [t' o t2) + (i, t' o hM)) ■ (4.2) 

We will get the required behaviour if there is a significant contribution to (4.2) 
from 1/3) = |0) and \a) ~ I7) for some \a). If this happens, the time-dependence 
of (4.2) is 

9{t - h)e{t' - t2) sin{(t - h){E^ - Eo)/h} sin{(<' - t2){E^ - Eo)/h}, (4.3) 

and then the Fourier transform has poles at the frequency 

noc. ^ {E^ - Eo)/h (4.4) 

which looks like what one gets from (3.1.7). But of course, as well as (4.3), 
there will be other contributions to (4.2), which are not functions of just the 
two variables (t — ti) and (t' — ^2)- So there is in general no reason to expect 
(4.4) to dominate. 

An heuristic quantum Langevin equation has been used for a system in- 
volving a shunted Josephson junction [3], [3], [TU], [T^], [H]. Then Q is the 
Josephson angle and 

ViQ) = -IQ-Iocos{Q) (4.5) 

with the constant coefhcients satisfying I > Iq. In this potential, there are no 
discrete quantum states, so the above considerations are irrelevant. What is 
more, in the classical motion Q increases indefinitely; but, when the damping 
due to the Langevin equation is included, the motion is in general qualitatively 
very different. So it is not possible to use an expansion in powers of the inter- 
action strength a (except perhaps for a short time interval after the initial time 
^o)- So in this case we are not able to study the validity of the approximations 
leading to the Langevin equation. 

Thus we are unable to find any conditions which would validate Senitzky's 
approximations for a general potential V{Q). 

5 Generalized Langevin equations 

Another route to a Langevin equation is via a generalized Langevin equation, 
constructed with the aid of projection operators, followed by the assumption 
that Q,P are "slow" variables (see for example [1], 0, [IS], [E], [H])- 

5.1 Projection operators 

In the dynamics defined in (2.1) and (2.2), we will assume that Q, P are 'slow" '(or 
'relevant') and qi,Pi are 'fast' variables. We need a projection operation, project- 
ing onto the 'slow' subspace. This entails a scalar product between operators, 
which we will choose to define as (for Hermitean operators A and B) 

(A, B)o = (1/2)((AB + BA))o = (l/2)tr{AB + BA)po}. (5.1.1) 
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This is symmetric which is convenient. (Some authors [5] use the more comph- 
cated Mori product. For simphcity, we do not take this course.) 

We have defined (5.1.1) using po, defined in (2.6). A disadvantage of (5.1.1) 
is that 

iA,[H,B])o^-{[H,AlB)o. (5.1.2) 

Another possibihty, which does not have this disadvantage, is to use the equi- 
hbrium distribution p. We will adopt this second choice in section 5.3 below, 
and then drop the suffices in (5.1.1) 

From now on, for simplicity, we will set to — 0. 

Given (5.1.1), the projection operation onto the "slow" subspace is defined 
to be, for any operator X, 

(X,Q(0))o JX^Pm^ 

- mmmo^^'^ + (p(o),p(o))o'^^'^' ^'-'-'^ 

and the complementary projection is (1 ^Vq). It seems to be usual to define Vq 
in terms of Q, P at the initial time, and we have emphasized this in (5.1.3). In 
general, a different choice of time would define a different projector. Note that 
(Q(0),P(0)) = because of invariance under time-reversal. 
Obviously 

^oQ'(O) ^ Q(0), (5.1.4) 

and for this reason the method may not be appropriate except in the linear case 
(2.3) (see section 8.2 of [2 ). Also 

VoQit)^Q{t) (t^O). (5.1.5) 

(For this reason, the projection method may be more appropriate in the Schrodinger 
picture, as in the derivation of the master equation.) But, for the linear case, 
we do have 

VoQoit) = Qoit), (5.1.6) 

and 

'Po9{qoi{t),Poi{t)) = 0, 'Po{QQ{t)g{qot{t),poi{t)} = Qo{t){g{qoi{t),PQi{t))Q 

(5.1.7) 

etc., where g is any function. These equations would not be true if we had 
defined V in terms of p rather than Vq in terms of pQ. 

5.2 Generalized Langevin Equation 

We need also to define the Liouville operator C, acting on an operator X, by 

CX^i[H,X]. (5.2.1) 

From the Heisenberg equations of motion, the generalized Langevin equation is 
deduced (see for example [2]): 

MQ{t) = P{t) (5.2.2) 
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P{t) + MQ'^Qit) - aAngQit) ~ aKKpP{t) 

+ / dt'{jQ{t - t')Q{t') + Jp{t - t')P{t')} = F{t), (5.2.3) 
Jo 

where 

F{t) = exp{i(l - Vq)C}{1 ~ Vq)CP{0) (5.2.4) 
with the property that 

VoF{t) = 0, (5.2.5) 

Kkq = (i^(0), Q(0))o/(Q(0), g(0))o, Kkp = {K{0), P{0))o/{P{0), P(0))o, 

(5.2.6) 

and 

jQ{t) = -(/:P(t),Q(0))o/(Q(0),Q(0))o, Jpit) = -(/:F(t),P(0))o/(P(0),P(0))o. 

.... ^^-^-^^ 

The quantities in (5.2.6) are in fact zero as we are using po in this section; but 
the corresponding quantities in section 5.3 (using p) are not both zero, 

One should note that (5.2.3) is an exact consequence of the Heisenberg equa- 
tions of motion, and yet it contains reference to the choice of po in (5.1.1). Any 
apparent dependence on po must cancel between the terms in (5.2.3). This, of 
course, may no longer be true if any approximations are made to (5.2.3). 

In (5.2.3), it is usual to think of F as being some generalized noise, and 
the J term as representing some (in general non-Markovian) friction. For this 
interpretation, we would hope that 

{F{t))o = 0. (5.2.8) 

We cannot prove this in general, or even that {F{t))o is time independent. (But 
in this connection, see section 8.3 of [2 .) 

It may throw some light on (5.2.3) to expand F and J through second order 
in a. Let 

F{t) ^ aFi{t) + a^F2{t) + .... (5.2.9) 

In (5.2.4), 

(1 - Vo)CP{0) = {l- Vo){~Mn^Q{0) + aK{0)} 

= (1 - Po)^o(O) + O(a') - Ko{0) + O(a'), (5.2.10) 
and so, from (5.2.4) and (5.2.7), to first order 

Flit) = exp{t(l - Vo)C}Ko{0). (5.2.11) 

Since we are working to first order in a here, we can neglect H' in C, and get 

(1 - Vo)CKo{0) ~ (1 - Po)i^o(O) = ii:o(0). (5.2.12) 

(using (5.1.7)), and similarly for higher derivatives. Therefore 

Fi{t) = exp(t/:)is:o(0) + O(a^) = Ko{t) + 0{a^). (5.2.13) 
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Inserting (5.2.13) into (5.2.7) we find a contribution (using (2.2)) 



jQ{t) ~ -{i/n)a{[H' ,Ko{t)],Q)o/{Q,Q)o - -{i/n)a^{[Ko{t),Kom}o, 

Jp{t) ~ 0, (5.2.14) 

which, working through second order, is all we need. 

We can now infer F2 in (5.2.9). The generalized Langevin equation (5.2.3) is 
exact, and our expansion through second order in a must be exact to that order. 
But the Langevin equation (3.1) with (3.6) is also exact through to second order 
(because the terms neglected in going from (3.1.3) to (3.1.4) were higher than 
the second). So (5.2.3) must be the same as (3.6) when (5.2.13) and (5.2.14) 
are inserted. This requires that 

F2{t) = [i/h) f dt'{[Ko{t),Ko{t')] - {[Ko{t),Ko{t')])o}Qo{t'). (5.2.15) 
Jo 

We may check that this satisfies (5.2.5). 

Thus we see that the 'noise' term (5.2.15) corrects the approximate friction 
term (5.2.14) so as to give the correct (through order a^) friction term in (3.6). 
The interpretation of F as noise may be open to question. 

5.3 The slow variable approximation 

It has been proposed to derive an ordinary Langevin Equation from the gener- 
alized one (5.2.3) as an approximation assuming that, in our example, Q and P 
are 'slow' variables. We will express this assumption in the form 

Q = 0{n)xQ, P = 0{n)xP (5.3.1) 

where Q is small. Presumably this means that <C w, where a) is some sort of 
typical value of the Ui in (2.2), 

We will follow the argument as presented in section 8.6 of [5]. First, we 
must depart from the choice po, which we have made up to now, and use the 
equilibrium distribution p as in [2j. Then, since [H,p] = 0, 

{X,CY) = -{CX,Y) (5.3.2) 

This allows (5.2.7) to be written in the form 

Jp(t) = (/(i)F(0),F(0))/(F(0),P(0)), 

jQ{t) = (l/M)(/(i)i^(0), (1 - P)P)/(g(0), Q(0)) = 0, 

where 

lit) = cxp{i(l - V)C}, 

and 

F(0) = (1 - V)CP{0) = (1 - V)i-n^Q{0) + aK{0)) 



(5.3.3) 
(5.3.4) 



16 



= a(l - V)K{0) = aK{0) - a''{cKQ/cQQ)Q{0) (5.3.5) 
where we define the functions of temperture 

CQQ = {Q{0),Qm, acKQ = iK{0),Q{0)). (5.3.6) 

There is now (using p) a nonzero contribution from (5.2.6) 

AxQ = acKQ/cQQ, Akp = 0. (5.3.7) 

the latter being a consequence of time-reversal invariance (and the symmetry of 
the scalar product). 

Note that here, because we have chosen the equilibrium distribution p and 
therefore can use (5.3.2), Jg = and Jp ^ 0; whereas in the approximation 
used in (5.2.13) (having chosen po) it was the other way round. 

For any X, (using (5.3.2)) 

^"^"^ - -(Q(o),Q(o))^^°) - {pwm ^'-'-'^ 

and if (5.3.1) is equivalent to 

CQ = 0{n) X Q, CP = 0{Q)xP (5.3.9) 

it follows that (5.3.8) is 0(fi); and then it is argued that we may approximate 
(5.3.4) by 

I{t) ~ exp{t£}. (5.3.10) 

In this approximation, 

F{t) ~ exp(i£)F(0) = aK{t) - {ckq / CQQ)Q{t) (5.3.11) 

and 

cppJp{t) ~ a'{K{t),K{Q)) - a\cKQ/cQQ){{K{t), Q(0)) + (Q(t), if (0))} 

+a\cKQ/cQQf{Q{t),Qm. (5.3.12) 

The last term in (5.3.11) (but not in (5.3.12)) is cancelled by (5.3.7). 

The approximations (5.3.11) and (5.3.12) are to be inserted into (5.2.2) to 
get the hoped for Langevin equation. We wil now make some remarks about 
this 'slow' approximation. 

(a) Because of the dependence on Q{t) on the right hand sides of (5.3.11) and 
(5.3.12), the resulting equation has a different form from an ordinary Langevin 
equation. 

(b) If we retain the term aK{t) on the right of (5.3.11) but discard everything 
else in (5.3.11) and (5.3.12), the resulting equation (5.2.3) would reduce to just 
the original Heisenberg equation of motion (3.1), not a Langevin equation at 
all. 
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(c) Since CP = — Jl^Q + aK, the assumption (5.3.9) seems to require that 
a as well as be small. 

(d) We can test the 'slow' approximation in the special case when X is a 
linear function 

= (5.3.13) 

when the ordinary Langevin equation is well known and is in fact exact. This 
equation amy be written in the form (5.2.3) with 

Jq = 0, Jp{t) = a'{l/M)Y,{c^l'm^uJ^fcos{uJ^t) = a' (MT) " Hr{/io (<)ifo (O)pis} , 

F{t) = aK^it) - a2Q(0) Y,{cjm,u:,f cos(c^,t), (5.3.14) 

where we have used the notation of (2.6). 

Let us see if we can get (5.3.14) from (5.3.11) and (5.3.12) in any approxi- 
mation. The nearest we can get is to write (5.3.11) and (5.3.12) as 

F{t) = aKo{T) + 0(a'^), 

Jp(t) = a2(l/Mr)(ifo(<)^o(0))o + 0(a3). (5.3.15) 

This does reproduce (5.3.14) except for the Q(0) term at the end of (5.3.14), 
which has to be considered as being of the neglected orders. 

This example confirms our expectation in (c) above, that a as well as Vl has 
to be regarded as small. But even then the 'slow' approximation seems to be 
incapable of reproducing the whole of the correct result (5.3.14) including the 
Q(0) term. 

It is worth remarking that the Langevin equation for linear model (5.3.13) 
can alternatively be written in the form (5.2.3) with Jp = and Jq ^ 0, 
and then there is no Q(0 term; but the slow approximation to the generalized 
Langevin equation forces the alternative form with Jq = 0. 

Finally, we note that if we are prepared to neglect higher powers of a on the 
right hand sides, we can immediately make Senitzky's approximation (3.5) for 
example, without going via the generalized Langevin equation. 

We conclude that the 'slow' approximation to the generalized Langevin equa- 
tion is not straightforward. 

6 Conclusions 

Although the Langevin equation (for a simple system interacting with a large 
environment) is often used, not so much attention has been given to judging its 
validity. We have critically examined two possible methods for doing this, with 
particular attention to the quantum case. 

The first is a direct approach due to Scnitzky [T]. He made approximations 
whose validity is not obvious in general. We identify one region in which the 
method may be justified, and that is where the measured frequency is close to 
the natural frequency of the free system. This condition cannot be formulated 
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when the system is nonhnear, although there are important examples where this 
is the case. 

The second approach (which also requires the system to be a linear one) 
seems to be completely different (see for example [5]). It proceeds by using a 
projection operator to a generalized Langevin equation, which is an exact con- 
sequence of the Heisenberg equations of motion, but which is of little practical 
use as it stands. In section 5.2, we argued that the identification of the noise 
and dissipation terms is ambiguous in the generalized equation. 

Then we studied the assumption that the system variables are 'slow' com- 
pared to the environment ones. This leads to an equation looking like an or- 
dinary Langevin equation. We examine the steps going into this derivation, 
particularly by comparison with the model in which all the equations of motion 
are linear, for which the Langevin equation is easily established and is exact. 
We argue that the 'slow' approximation necessarily entails neglect also of terms 
of higher order in the coupling strength a, as well as in the frequency ratio. But 
even then the correct Langevin equation requires a selective choice of orders of 
a. 

JF would like to thank CNPq, Brazil, for a grant. 
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